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Abstract. We explore the border between regular and chaotic quantum dynamics, 
characterized by a power law decrease in the overlap between a state evolved under 
chaotic dynamics and the same state evolved under a slightly perturbed dynamics. 
This region corresponds to the edge of chaos for the classical map from which the 
quantum chaotic dynamics is derived and can be characterized via nonextensive entropy 
concepts. 



1 Introduction 

Classical chaotic dynamics is identified by extreme sensitivity to initial condi- 
tions. Under chaotic dynamics, two arbitrarily close points in phase space diverge 
at an exponential rate, quantified by the Lyapunov exponent [1]. Non-chaotic 
dynamics does not show this extreme dependence on initial conditions, hence, 
the Lyapunov exponent is equal to zero. However, it has been conjectured, that, 
though the Lyapunov exponent may vanish, the system may have a positive gen- 
eralized Lyapunov coefficient [2] describing power-law, rather than exponential, 
divergence of classical trajectories. This is the case, at the border between chaotic 
and non-chaotic dynamics (the 'edge of chaos') where the Lyapunov exponent 
goes to zero, but there remains a positive generalized Lyapunov coefficient. 

This paper identifies a characteristic signature for the edge of quantum chaos. 
Quantum states maintain a constant overlap fidelity (heretofore referred to as 
the overlap or fidelity) , or distance, under all quantum dynamics, regular and 
chaotic. One way to characterize quantum chaos is to compare the evolution of 
an initially chosen state under the chaotic dynamics with the same state evolved 
under a perturbed dynamics [3-5]. When the initial state is in a regular region 
of a mixed phase space system, a system whose phase space has regular and 
chaotic regions, the overlap remains close to one. When the initial state is in a 
chaotic region, the overlap decay is exponential. This paper explores the edge of 
quantum chaos, a region of polynomial overlap decay [6]. 
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This paper is structured as follows, we first give a short review of the Lya- 
punov description of chaotic dynamics and the lack of correspondence between 
this description and quantum dynamics. We then discuss suggested character- 
istics of quantum chaos (known as signatures of quantum chaos), concentrating 
on the overlap decay first introduced by Peres [3,4] . Overlap decay proves to be 
a useful signature of quantum chaos from which to explore the border between 
regular and chaotic quantum dynamics. Next, we will review the nonextensive 
entropy form introduced in [7] and discuss its relevance to the edge of chaos 
phenomenon. Finally, we locate the 'edge of quantum chaos' in a mixed phase 
space system and, using the nonextensive entropy formalism, show how the over- 
lap decay at the edge of quantum chaos depends on perturbation strength and 
Hilbert space dimension. 



2 Classical Chaos 

The Lyapunov exponent description of classical chaos is as follows [1]. Let Ax be 
the distance between two points on phase space. We define £ = faAi(o)^o( ^(o) )> 
to describe how far apart two initially arbitrarily close points on phase space be- 
come at some time t. Generally, is the solution to the differential equation 

giving the solution 

m = e Al * (2) 

where Ai is the Lyapunov exponent (the use of the subscript will become clear 
later on) . As seen from the above equation when the Lyapunov exponent is pos- 
itive two arbitrarily close points on phase space diverge at an exponential rate. 
Thus, the dynamics described by £(i) is strongly sensitive to initial conditions 
and we have chaotic dynamics. 



3 Quantum Chaos 

While the Lyapunov exponent description of chaos works well for points on 
a classical phase space it cannot hold true for quantum mechanical states or 
wavefunctions. A measure of distance between quantum wavefunctions is the 
overlap 

O t = (#|«P>. (3) 

However, the overlap between two quantum wavefunctions remains unchanged 
under unitary evolution governed by the linear Schrodinger equation. This is 
seen from 

0(n) = (*{U n yU n \$) = Oi, (4) 

where U is the unitary system evolution. Hence, the distance between two arbi- 
trarily close quantum mechanical wavefunctions, like the distance between two 
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Liouvillc probability densities, docs not diverge and cannot be described by the 
Lyapunov exponent picture. This seeming lack of correspondence has led to the 
study of 'quantum chaos,' the search for characteristics of quantum dynamics 
that manifest themselves as chaotic in the classical realm [8-11]. 

Many such characteristics, quantum signatures of chaos, have been detailed 
in the literature and tested for quantum analogs of classically chaotic systems. 
These signatures can be divided into two broad categories, static signatures and 
dynamic signatures. Static signatures look at characteristics of the Hamiltonian 
or unitary operator governing the system. The conjecture is that the evolution 
operator of quantum chaotic systems have statistical properties similar to those 
of random matrices. Hence, quantum analogs of classically chaotic systems show 
level repulsion, that is, if the energy eigenvalues of the system are ordered, 
the difference between nearest neighbors would result in a histogram with a 
Wigner-Dyson distribution [11] and not of a Poisonnian distribution expected for 
regular systems. In addition, the magnitude of the elements of the eigenvectors of 
quantum chaotic operators follow xl distributions[12] from appropriate random 
matrix ensemble. 

Dynamic signatures of quantum chaos look at the evolution of a state un- 
der the quantum chaotic operator compared to the same evolution with some 
additional perturbation. An example of a dynamic signature of chaos is hyper- 
sensitivity to perturbation. For chaotic systems (both classical and quantum), 
the amount of information necessary to track the state of a system when the 
dynamics is interrupted by an unknown perturbation grows at an exponential 
rate with increasing time [5]. Other signatures of quantum chaos look at the 
entropy growth of chaotic systems versus regular systems [13]. 

4 Overlap Decay 

Overlap decay as a signature of quantum chaos was first introduced by Peres 
[3,4,14] as a quantum analog of initial state sensitivity. Rather than look at 
two slightly different states and see how they evolve under a certain dynamics, 
Peres suggested looking at one state, |V>i), and see how it evolves under under 
two slightly different dynamics, an unperturbed Hamiltonian H , and the same 
Hamiltonian with a small perturbation H + 5V, where S is the perturbation 
strength. The overlap at time t is 

0{t) = \(Mt)\Mt))\- (5) 

where \tp u (t)) = eT %m \i\>i) and \i> p (t)} — e^ l( - H+5v ^ t \ipi} are the unperturbed and 
perturbed states, respectively. The initial behavior of the overlap shows different 
behavior depending on whether or not H is chaotic. Recently, the study of overlap 
decay has seen a revival of interest which has served to identify several regimes 
of overlap decay behavior based on whether the system is chaotic or regular, the 
type and strength of perturbation, and the type of initial state. We note that 
many of the works cited use O 2 as the fidelity. Here, we follow [15] and simply 
use the overlap, O. 



4 



Yaakov S. Weinstein et al. 



For quantum chaotic systems, quantum versions of classically chaotic maps 
or random matrix models, there are several regions of behavior based on the 
strength of the perturbation. Perturbation strength is characterized by the size 
of a typical off diagonal clement, a, of the perturbation operator when the per- 
turbation operator is written in the ordered eigenbasis of the unperturbed dy- 
namics, H. If a is less than the average level spacing of the unperturbed system, 
A, the perturbation is weak. If a > A, the perturbation is in the Fermi Golden 
Rule (FGR) regime [16,17]. The average level spacing is equal to 2-k/N where N 
is the dimension of the system Hilbert space. A typical off diagonal element of 
the perturbation operator in the ordered eigenbasis of the unperturbed system 
(where the unperturbed system is assumed to have eigenvector statistics of a 



random matrix see [18]) is equal to y <5 2 V^„, where V^ n is the second moment 
of the matrix elements of the perturbation Hamiltonian. 

For short enough time, the overlap decay is quadratic for any perturbation 
strength [3]. After this time, weak perturbations lead to a Gaussian overlap 
decay as expected from first order perturbation theory [14]. Perturbations in the 
FGR regime lead to an exponential decay of overlap whose rate, r, increases 
with perturbation strength. For many systems, r increases as the perturbation 
strength squared [16], however, some systems do not show this exact dependence 
[19]. For systems with a classical analog, initial coherent states, and classical 
perturbations, r will increase with perturbation strength until the decay rate 
reaches a value given by the Lyapunov exponent of the analog classical system 
[16,20,21]. This is a wonderful example of a dynamical property of the classical 
system emerging in quantum mechanics. The increase r may also saturate at 
the bandwidth of H [16]. 

For initial states that are coherent or random states, the Gaussian or expo- 
nential behavior continues until the fidelity reaches a saturation point at ~ 1/N. 
For initial states that are system eigenstates, there is no decay in the weak per- 
turbation regime, and in the FGR regime there is an exponential decay of overlap 
which saturates at ~ 1/T [19,22]. 

Regular, non-chaotic, systems have a Gaussian decay in the FGR regime [15]. 
The Gaussian decay is faster then the exponential decay of chaotic states. This 
result has been explained using correlation functions and may be understood as 
follows: a perturbation to a chaotic system is quickly spread out to the entire 
Hilbert space of the system. Repetition of the same perturbation does not lead to 
a compounded error (the correlation time is short). Regular systems, however, do 
not have this mechanism to spread out the perturbation. The same perturbation 
compounds the error (there is a long correlation time) leading to a faster fidelity 
decay. For some more complex regular systems, a power-law decay has also been 
observed [23]. 

Here we study a mixed system, a system with both chaotic and regular 
regimes. Coherent states within the regular regime are practically eigenstates 
of the system and the overlap of these states oscillates close to unity [4,15], as 
shown in figure 1. Coherent states in the chaotic regime behave like other chaotic 
systems, they show exponential overlap decay in the FGR regime and Gaussian 
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overlap decay for weak perturbations. We explore the border between the regu- 
lar and chaotic regions and show that, in both the FGR and weak perturbation 
regimes, coherent states near this border have a polynomial overlap decay [6]. 




200 400 600 t 800 



Fig. 1. Overlap decay of coherent states in the regular and chaotic zone of the J = 480 
QKT with k = 3 and 8 = .005 (in the FGR regime) . The overlap of the coherent state 
centered at the fixed point of the classical kicked top is practically an eigenstate of the 
QKT and oscillates close to unity. The overlap of the coherent state in the chaotic zone 
decays first quadratically and then exponentially. 



5 Nonextensive Entropy 

The Boltzmann-Gibbs formulation of statistical mechanics is one of the pillars of 
modern day physics. However, its applicability indeed its formulation, is only for 
systems with short range interactions. Self-gravitating systems and other systems 
such as those with long range forces, long range memories (non-Markovian) , or 
multifractal structures present significant difficulties for the application of the 
Boltzmann-Gibbs formulation. In 1988, in an attempt to understand the nature 
of some of these systems, one of us [7] proposed the generalization of Boltzmann- 
Gibbs statistical mechanics on the basis of a generalized non-extensive entropy, 
S q . When q = 1, S q recovers the usual Boltzmann-Gibbs entropy. 
The generalized entropy is given as follows: 

1 - T W v q 

S q = k ^ (6) 

where k is a positive constant, pi is the probability of finding the system in 
microscopic state i, and W is the number of possible microscopic states of the 
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system; q is the entropic index which characterizes the degree of the system 
nonextensivity. In the limit q — > 1 we recover the usual Boltzmann entropy 

w 

Si = -k^pi Inpi. (7) 

i=l 

To demonstrate how q characterizes the degree of nonextensivity of a system, 
we present the S q entropy addition rule [24]. If A and B are two independent 
systems such that the probability p( A + B) = p(A)p(B), the entropy of the total 
system S q (A + B) is given by the following: 

SM + B) S q (A) S q (B) S q (A)S q (B) 

k -^- + ^- + (1 - q > P • (8) 

From the above equation one can see that q < 1 corresponds to superextensivity, 
the combined system entropy is more then the sum of the two independent 
systems; q > 1 corresponds to subextensivity, the combined system entropy is 
less then the two independent systems. Using this entropy to generalize statistical 
mechanics and thermodynamics has helped to explain many natural phenomena 
in a wide range of fields. 

One application of nonextensive entropy occurs in one dimensional dynam- 
ical maps. As explained above, when the Lyapunov exponent of a system is 
positive, the system dynamics is strongly sensitive to initial conditions and 
is called chaotic. When the Lyapunov exponent goes to zero it has been con- 
jectured [2] (and proven for the logistic map [25]) that the distance between 
two initially arbitrarily close points can be described by ^ = A 9s(!n £ 9se ™ lead- 
ing to £ = [1 + (1 — qsen)^q scn t] 1 ^ 1 ~ qs "™*> (sen stands for sensitivity). This 
requires the introduction of X q3en as a generalized Lyapunov coefficient. The 
Lyapunov coefficient scales inversely with time as a power law instead of the 
characteristic exponential of a Lyapunov exponent. Thus, there exists a regime, 
Qsen < 1, Ai = 0, Aq 3e7i > 0, which is weakly sensitive to initial conditions and is 
characterized by having power law, instead of exponential, mixing. This regime 
is called the edge of chaos. 

The polynomial overlap decay found for initial states of a mixed system 
near the chaotic border are at the 'edge of quantum chaos', the border between 
regular and chaotic quantum dynamics. This region is the quantum analog of 
the classical region characterized by the generalized Lyapunov coefficient. 

6 The Quantum Kicked Top 

The system studied in this work is the quantum kicked top (QKT) [26] defined 
by the operator: 

Uqkt = e - inJ y/ 2h e- iaJ */ 2Jh . (9) 

J is the angular momentum of the top and a is the 'kick' strength. The rep- 
resentation is such that J z is diagonal. The classical version of the kicked top 
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has either regular, mixed, or fully chaotic dynamics depending upon the kick 
strength. We work with a QKT of a = 3 whose classical analog has a mixed 
phase space, with clearly defined regions of chaotic and regular dynamics. The 
perturbed operator is simply a QKT with a stronger kick strength a'. Hence, 
the perturbation operator, V = Sir J 2 /2 J where 8 = a! — a. 

When J is even the QKT has three invariant, dynamically independent sub- 
spaces [4,26]. We can write basis functions for the invariant subspaces in terms 
of the eigenvectors of J z which, following Peres [4], we will denote as \m). The ee 
subspace is even under a 180° rotation about x and even under a 180° rotation 
about y and has a Hilbert space dimension of TV = J/2 + 1. The basis functions of 
the ee subspace are |0) and (|2m) + | — 2m))/-\/2, where m ranges from 1 to J/2. 
The oo subspace is even under a 180° rotation about x and odd under a 180° 
rotation about y and N = J/2. Its basis functions are (|2m — 1) — |1 — 2m))/v / 2- 
The oe subspace is odd under a 180° rotation about x, N — J, and has basis 
functions (|2m) - 1 -2m))/V2 and (|2to-1) + |1-2to))/\/2. All of the numerical 
simulations in this work were done in the oo subspace of the QKT. This is to 
say, we construct the complete QKT and transform it into a basis such that it 
is block diagonal with the dimensions of the three blocks mentioned above. The 
columns of the transformation matrix T to get the QKT operator into block 
diagonal form are the states which form the bases of the invariant subspaces. 
We take only the block corresponding to the oo subspace and use that as our 
map. The initial angular momentum coherent states are also transformed into 
this basis and, again, only the part of the state corresponding to the oo subspace 
is used. 

The phase space of the classical kicked top is the unit sphere, x 2 + y 2 + z 2 = 1 
and the resulting action of the map is: 

x' = z 

y' = x sin(az) + y cos(az) (10) 
z' = —x cos(az) + y sin(az). 

For a = 3 there are two fixed points of order one at the center of the regular 
regions. They are located at 

Xf = Zf = ±0.6294126, y f = 0.4557187. (11) 

The regular and chaotic regions of the kicked top are clearly seen in the classical 
maps phase space shown in the figure 2. Hence, we expect quantum coherent 
states centered near the classical periodic point to exhibit significantly different 
behavior than coherent states centered in the classically chaotic region of the 
map. This is indeed seen in figure 1. 

7 Locating the Edge of Quantum Chaos 

To locate the edge of quantum chaos we use initial angular momentum coherent 
states keeping y equal to yf of the positive fixed point for the classical kicked 
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top and changing z until a state which has a power-law overlap decay is found. 
The state which gives this behavior depends on the angular momentum of the 
QKT, J, but, given a fixed J, the power law emerges for perturbations in both 
the weak perturbation and FGR regimes. Examples of edge of quantum chaos 
behavior in both regimes is shown in figure 3 for a QKT of J = 240. As seen 
in the figure, the power-law overlap decay is transitory between the quadratic 
and exponential behavior of the overlap decay. This transitory region does not 
appear for chaotic states (as shown in figure 1) or states close to the fixed point 
of the map, it is a signature of the 'edge of quantum chaos.' 

The overlap decay for the edge of quantum chaos state is very well fit by the 
solution of the differential equation 

dO/d(t 2 ) = -0^/ T l er (12) 

In the above equation, rel stands for relaxation which is an appropriate descrip- 
tion of a q > 1 phenomenon. In classical systems, q re j characterizes the relaxation 
of initial states towards an attractor. Although we do not know how to derive 
this differential equation from first principles, the numerical agreement is re- 
markable (see also [27]). A timc-dependent g-cxponcntial expression analogous 
to the one shown here has recently been proved for the edge of chaos and other 
critical points of the classical logistic map [25] . 

The values of q re i and r qrel depend on the perturbation strength, J, and on 
whether the perturbation strength is in the weak perturbation or in the FGR 
regime. For perturbation strengths below the FGR regime, 5 is less than some 
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Fig. 2. 10000 points of a chaotic orbit on the phase space of the classical kicked top 
with kick strength a = 3. The regular regions are clearly visible. The spherical phase 
space is projected onto the x — z plan e by multi plying the x and z coordinates of each 
point by R/r where R = ^2(1 - \y\) and r = ~y 2 ) [4]. 
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Fig. 3. Overlap versus time for an initial angular momentum coherent state located 
at the border between regular and chaotic zones of the QKT of spin 240 and a = 3. 
This region, the edge of quantum chaos, shows the expected power law decrease in 
overlap. The top figure is for a perturbation strength in the weak perturbation regime, 
S — .0003 and the bottom figure is for a perturbation strength of S = .01, within 
the FGR regime. On the log-log plot the power law decay region, from about 600- 
2500 in the weak perturbation regime and 20-70 in the FGR regime, is linear. We 
can fit the decrease in overlap with the expression [1 + (g re ; — l)(t/r qrcl J 2 ] 1 /^-?^) 
where, in the weak perturbation regime, the entropic index q re i = 3.3 and T Qrel = 1300 
and in the FGR regime q re i — 4.25 and T Qrel = 34. The insets of both figures show 
ln qrel O = (0 1 ~ qr "' — 1)/(1 — q r ei) versus t 2 ; since ln 9 x is the inverse function of 
e q = [1 + (1 — q) x] !-« , this produces a straight line with a slope — 1/r 2 (also plotted). 



S c , q re i remains constant at a value of q c rel . The transition into the FGR regime 
arises when the typical off diagonal elements of V are larger than A. We can 
approximate S c ~ [16] where, in the oo subspace of the kicked top 

N = J/2. Values for 5 C are shown in the above table. In the FGR regime q re i 
continues to increase with perturbation strength until a saturation perturbation, 
S s , after which q re i — q^. el remains constant. As the top becomes more classical, 
with increased J, there is a decrease in q£ el while in the weak perturbation 
regime, but an increase in the rate of increase of q re i while in the FGR regime. 
Because of this larger rate of increased q re i in the FGR regime, q^. el increases with 
increasing J. However, we see that this is not the case for the J = 120 case. For 
J = 120, q re i increases beyond the expected saturation point. This may be due 
because the J = 120 coherent state is so large that, at stronger perturbations, 
it leaks out of the regular region of the map in more than one place due to the 
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Table 1. The edge of quantum chaos, the critical perturbation strength, S c , and q^ el , 
the q re i for perturbations strengths below 5 C , for explored values of J. As J increases 
behavior characteristic of the edge of quantum chaos occurs further away from the fixed 
point of the classical map. The critical perturbation and q c rel decrease with increased 
J. 
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odd shape of the regular region (see figure 6). The value of r qrel decreases with 
perturbation strength and is well fit by a line of slope approximately -1 on a log- 
log plot; T qrcl also decreases with increasing J at a fixed perturbation strength. 
The values of q re i and T qrcl for a number of different perturbation strengths can 
be seen in the figure 4. 

The location of states exhibiting edge of quantum chaos behavior is not the 
same as the edge of chaos for the classical kicked top. This is due to the finite 
size of the coherent states. Though the coherent state may be centered at a point 
that is classically regular, some of the state may 'leak out' into the chaotic region 
of the map. As J increases the coherent state gets smaller and the quantum edge 
approaches the classical one. The location of the edge and the size of the coherent 
state compared to the regular regions of the map are shown in figure 6. 

In conclusion, we have explored the region at the between chaotic and non- 
chaotic quantum dynamics, the edge of quantum chaos. Coherent states located 
at this border exhibit a power-law decrease in overlap as opposed to practically 
no decay for coherent states near the periodic point of the classical map and 
exponential overlap decay exhibited by fully chaotic quantum dynamics. This 
region is the quantum parallel of the classical region at the border between 
regular and chaotic classical dynamics where the Lyapunov exponent goes to 
zero and the mixing is characterized by the generalized Lyapunov coefficient. 
Further studies of this rich system are certainly welcome. 

One of us (C.T.) acknowledges warm hospitality by H.-T. Elze and the or- 
ganizers during the interesting meeting. 
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Fig. 4. Values of q re i and r q (inset) for J = 120 (x), 240 (circles), 360 (diamonds) 
and 480 (stars). q re i remains constant for perturbation strengths below the critical 
perturbation and above the saturation perturbation. In between q re i increases with a 
rate dependent on J. The values of q^ei, Qreh <^c and S 3 can be seen in the figure. In 
addition the rate of growth of q re i with increased perturbation strength can be seen. 
The inset shows a loglog plot of the value of r q versus 8 for the above values of J. The 
data can be fit with a lines of slope -1.06, -1.03, -1.07, and -1.08 for J = 120, 240, 360 
and 480 (top to bottom). 

4i 1 1 1 1 1 

3.8 

3.6 

3,4 

D" 

3.2 
3 
2.8 



.002 .004 .006 .008 1/J 



Fig. 5. Values of q^ el versus 1/J. These are determined by exploring a number of 
perturbations much less than 8 C . 
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Fig. 6. Coherent states of J = 120 (circles) and J = 480 (x) at the edge of quantum 
chaos superimposed on the upper right quarter of the classical phase space (shown in 
figure 2). The spherical phase space and ellipsoidal coherent states are projected onto 
the x — z plane as described above. The J = 120 coherent state is significantly larger 
then the J = 480 coherent state, causing edge of quantum chaos behavior to appear 
for the J = 120 wavefunction much closer to the classical periodic point than edge of 
chaos behavior for the J = 480 wavefunction. The inset is a semilog (y-axis) plot of 
1/J versus distance from the fixed point of the classical map. Note that the point of 
appearance of edge of chaos behavior follows an inverse power law with the size of the 
coherent state. 
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